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Starting from Plebanski formulation of gravity as a constrained BF theory we 
propose a new spin foam model for 4d Riemannian quantum gravity that generalises 
the well-known Barrett-Crane model and resolves the inherent to it ultra-locality 
problem. The BF formulation of 4d gravity possesses two sectors: gravitational and 
topological ones. The model presented here is shown to give a quantization of the 
gravitational sector, and is dual to the recently proposed spin foam model of Engle 
et al. which, we show, corresponds to the topological sector. Our methods allow us 
to introduce the Immirzi parameter into the framework of spin foam quantisation. 
We generalize some of our considerations to the Lorentzian setting and obtain a new 
spin foam model in that context as well. 



I. INTRODUCTION 

The spin foam approach to quantum gravity was born about a decade ago as a result of 
cross-fertilization between several different fields: loop quantum gravity, topological quan- 
tum field theory and simplicial gravity. On one hand, this approach is supposed to give 
a spacetime covariant version of the picture of quantum geometry arising in loop quan- 
tum gravity. On the other hand, it attempts to apply methods and ideas that are known 
to be very powerful in the context of 3-dimensional quantum gravity to the problem of 
quantization of general relativity (GR) in 4 spacetime dimensions. At the same time, the 
models considered are in the simplicial setting, and are either derived from or inspired by a 
discretization of some classical gravity action. 

The most well-known model -that due to Ponzano and Regge [H- appeared in the sixties 
in the context of three-dimensional quantum gravity, thirty years before the term "spin foam 
models" and twenty years before the first works on loop quantum gravity. In the context of 
four spacetime dimensions the notion of (what we now call) a spin foam as a history of spin 
network evolution was introduced, from the perspective of loop quantum gravity, in 0,11]. 
A new era for the subject was opened by an important work by Barbieri who realized 
that quantum states of loop quantum gravity can be understood (and heuristically derived) 
by applying the quantization procedure to the usual geometric tetrahedron in 3 spatial 
dimensions. Even though this work may have not received a proper credit in the literature, 
it was extremely influential in that it suggested that at least to some extent quantum gravity 
may be about quantizing geometric structures. This idea was then quickly applied to the 
problem of quantization of a geometric simplex in M^, with a breakthrough result being the 
construction of the now famous Barrett-Crane model [5|. A relation between this model and 
various versions of the simplicial gravity was explored in works [6| , 0] and [sl showing that 
the Barret-Crane model is related to a formulation of gravity as a constrained BF theory. 

The key notion arising in the Barrett-Crane model is that of simple (or class I) represen- 
tations of the Lorentz group in four dimensions. In work by Freidel, Krasnov and Puzio 
these representations were understood as coming from the harmonic analysis on the certain 
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homogeneous group space, a crucial fact that opened a way to the Lorentzian generahzation 
Tol of the original Barrett-Crane model in Euclidean signature. 



As another byproduct to these developments came a realization that spin foam models 
of the Barrett-Crane type are naturally obtainable from a quantum field theory on a certain 
group manifold This generalizes the field theory introduced by Boulatov [l2[ in the 

context of 3-dimensional quantum gravity. This way of deriving spin foam models later 
received a catchy name of "group field theory" , see [131] for a review. These developments 
were important at least for the fact that they allowed to have an efficient bookkeeping tool 
for the various amplitudes that are part of the spin foam model (the amplitudes for the 
lower-dimensional "cells" were not very clearly understood and/or described in the work 
of Barrett and Crane j^). Most optimistically, the group field theory method may have 
a fundamental status and has a potential to lead to an independent approach to quantum 
gravity. 

After this initial period of rapid development the subject of spin foam models entered a 
much calmer stage. During this stage evidence was accumulated that casted serious doubts 
on the physical correctness of the Barrett-Crane model(s). Namely, it became clear that, 
in spite of correctly capturing at least some aspects of quantum geometry of a single sim- 
plex in 4 dimensions, this theory is likely not viable as a theory of quantum geometry of 
more complicated objects constructed from several such simplices. The suggested problems 
with the models were of two main types: 1) the absence or very limited nature of correla- 
tion between amplitudes of neighboring simplices glued together (the so-called ultra-locality 
problem); 2) analytical and numerical evidence (in the case of a single simplex) that the 
simplicial geometries dominating the partition function of the model in the limit of large 



representations are the degenerate ones [IJ], [15|, [16|. These problems make it unlikely 



that the usual world described by a non-degenerate metric can arise in the semi-classical 
limit of these models. Many other criticisms (as well as some counterarguments to them) 
of the Barrett-Crane model(s) were brought forward, and for the past years the subject 
experienced some stagnation. 

A new development came with the work by Livine and Speciale [13], who used the coherent 
state methods to rewrite the formulae for the BF theory vertex amplitudes. They also 
suggested that these methods may prove instrumental in finding an improved version of the 
Barrett-Crane model that is free from the above mentioned problems. 

Shortly afterwards, there appeared an important paper by Engle, Pereira and Rovelli [l8| , 
who proposed a new model that was by construction free from the ultra-locality problem. 
The authors of [l8| argue that the key problem of the Barrett-Crane model(s) is that they 
impose certain geometric constraints (the so-called cross-simplicity constraints) "strongly" , 
while some other approaches to the problem, notably the canonical approach, suggest that 
these constraints should be imposed only weakly, if one is not to loose the degrees of freedom 
of the gravity theory in question. On a more technical side, the paper [l8| also hinted at a 
new way in which the above mentioned geometric constraints can be looked at. This new 
interpretation of the "simplicity" constraints was further developed by Alexandrov in his 



recent wor k 1 191 ]. While the present work was in its final stages, an extended version [20| of 
the paper 3] has appeared. 

The present work resulted from an attempt to understand and extend the construction of 



181 ] using methods similar to [17|. To this end, we have further developed the new geomet- 



ric perspective on the simplicity constraints, as well as the new geometrical interpretation 



suggested by the use of coherent states in [17[ . As the first immediate result of our analysis 
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came the realization that the problem of ultra-locality common to all Barrett-Crane-type 
models results from a way these models impose the so-called simplicity constraints. This 
way of dealing with the constraints focuses on a single 4-simplex and completely ignores the 
fact that the amplitude associated to this 4-simplex is to become a part of a quantum am- 
plitude associated to the whole triangulation. The use of coherent state makes this feature 
of Barrett-Crane models crystal clear, as well as suggests a way to remedy this problem. 

As a further result of our analysis came a realization that there is not one, but two distinct 
ways to impose the simplicity constraints. The fact that there should be two different 
distinct approaches to "simplicity" (as suggested by the fact that there are two distinct 
symplectic structures on the phase space in question) was noted long ago, and was one of 
the early arguments brought forward against the Barrett-Crane models, which left no place 
for a second description. This ambiguity gave grounds to suspect that the Barrett-Crane 
models might have given a quantization of the "wrong" sector of the theory. However, later 
some rather convincing arguments were put forward to the opposite effect, see in particular 



21[ . The coherent-state-inspired way of imposing the simplicity constraints leads to a much 
clearer understanding of this issue. Thus, we will see that there is not one, but two distinct 
models that should be interpreted as giving quantization of the gravitational and topological 
sectors of the theory respectively. One of the models that results from our analysis is that 
proposed in [l8|. The other model that we find is new. It is similar to the model of [l8| in 
the sense that the simplex amplitudes are strongly correlated and there is no ultra-locality. 
However, it also captures some of the main positive features of the Barrett-Crane model. The 
simplex amplitudes for this model are quite different from the model of Engle et al. [3] , see 
below for a detailed comparison. Importantly, we show that the model corresponds to 
the topological sector of the theory, while the new model we propose gives the gravitational 
sector. Our arguments are based on a detailed study of the gravitational and topological 
sectors and on an observation that the new model keeps certain key characteristics of the 
Barrett-Crane model. The model we propose thus results by merging some key features of 
the Barrett-Crane model with the new insight about ultra-locality achieved in fl8| . 

The model we obtain (as well as the model of 18j) has two main advantages as compared 
to the Barrett-Crane one: First, it clearly possesses a higher degree of correlations among 
the neighbouring simplices, thus overcoming the ultra-locality of the Barrett-Crane model; 
Second it uses SU(2) structures that bear more resemblance to those appearing in traditional 
loop quantum gravity. Moreover, the fact that we now have two distinct models that clearly 
correspond to two different sectors of the theory - the gravitational and the topological ones 
- removes the above mentioned shortcoming of the model of Barrett and Crane that failed 
to describe the opposite sector. This availability of the description of both sectors is what 
makes the new developments particularly exciting. 

Last but not least, one of the main outcomes of our analysis in the present work is that 
the coherent state methods suggest a way of naturally incorporating the so-called Immirzi 
parameter into the framework of spin foam models, something that was an open problem for 
a long time. Indeed, the absence of a Barrett-Crane-like model that would incorporate the 
Immirzi parameter was long thought to be a suspicious feature of the spin foam framework. 

We also describe very briefly how the new models can be obtained within the framework 
of group field theory. In addition, our construction allows us to generalize the arguments 
to the Lorentzian case and obtain a new model in this setting. The new Lorentzian model 
that we obtain has similar advantages as compared to [lo| . 

The organization of the paper is as follows. Section [TTl develops a new geometric perspec- 
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tive on the simplicity constraints. Section IIIII introduces the central technique of coherent 
states and gives the description of the topological BF theory model using this technique. 
We re-derive the Euclidean models of and [l8| in section HV] and also obtain a new model 
here. Section |V] describes how the Immirzi parameter can be introduced into the spin foam 
model framework. We generalize the constructions to the Lorentzian setting in section IVII 
and derive a new Lorentzian model. Section IVlII discusses how the models considered in this 
paper can be casted into the group field theory framework. In section IVIIII we discuss an 
alternative way to impose the constraints. We finish with conclusions. 

After a first version of this pap er was placed on the arxiv, a whole set of related works 
appeared. In particular, paper 22] is very close in spirit to the present work and paper 2^ 
studies further the finite Immirzi parameter case. 



II. A NEW APPROACH TO THE GEOMETRIC SIMPLICITY CONSTRAINTS 

Historically, the spin foam models were obtained by applying the quantization procedure 
to a geometric simplex in 4 dimensions. One starts by associating a Lie-algebra valued 
variable Xf to each of the 10 faces / of the simplex. The Lie-algebra in question is that of 
the Lorentz group appropriate to the situation at hand, thus either SO (4) in the Euclidean 
case or SL(2, C) in the Lorentzian case. A norm of the Lie-algebra element introduced this 
way receives the geometric interpretation of the area of the corresponding face (note that, 
in fact, there are two different norms on the Lie- algebras in question; the duality exchanging 
these two norms will play an important role in what follows). The quantization procedure 
then promotes this Lie-algebra element into a quantum operator, see more on this below. 
The spin foam models are then obtained by analyzing what certain geometrical constraints 
appropriate for a flat geometric 4-simplex boil down to at the level of the operators. Selecting 
the states that satisfy the operator version of these constraints one derives the simplex 
amplitude. 

It turns out to be possible, however, to get the same models by deriving them from a 
certain simplicial action that can be thought of as a precise discrete analog of the gravity 
one. The main idea for this comes from the so-called BF formulation of gravity theories, 
and from the systematic procedure of deriving spin foam models from the action principle 
proposed in [1]. 



A. BF formulation of gravity and path integral quantization 

In Plebanski (or BF, as it was subsequently called) formulation of gravitational theory 
3, H, 0, one starts by the ^0(4) (or 50(3, 1) in the Lorentzian case) BF theory 



J B^-^F{A)u (1) 



where F{A) is the curvature of an 5*0(4) connection and S is a Lie-algebra valued two-form 
field. This two form needs to satisfy certain "metricity" (or "simplicity") constraints 



(2) 

that guarantee that it is of the form either B^"^ = ±6^ A 6''^ (called sector I in 0) or 
B^-^ = ±e'-^^^^9K A Ol (called sector II in 0). Here 9^ are the tetrads, I,J,K,L are 
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"internal" indices taking values 0, 1, 2, 3, and e^-^^^ is the completely anti-symmetric tensor. 
Sector II gives the Cartan-Weyl form of Einstein's theory and we will refer to it as the 
gravitational sector. Sector I gives a theory with no local degrees of freedom. We will call 
this sector topological. 

The path integral quantization of gravity in the BF formulation has been studied in 
detail in [2^. The analysis is complicated by the fact that the dynamical system is question 
is quite non-trivial. Indeed, to understand the correct gauge-fixing terms that need to be 
added to the action to make sense of the path integral it is necessary to understand the 
structure of the symmetry algebra of the theory, and thus understand the structure of the 
phase space. Here, on top of the usual Gauss and the above simplicity constraints as primary 
constraints one finds also many secondary constraints (in the topological sector of the theory 
there are also tertiary constraints). The resulting constraint algebra has both first-class and 
second-class parts. The first-class constraints are to be gauge-fixed, which converts them 
into second class. The path integral can then be taken over the reduced phase space. It is 
important to take into account the fact that the path integral measure includes a square 
root of the determinant of the Dirac symplectic structure on the reduced phase space, and 
thus may be non-trivial. This non-trivial procedure has been carried out in Section IV of 
26| . and the resulting measure on the original space of fields has been found. 

However, what is often not appreciated is that the result of this complicated procedure is 
rather simple. If one, as we are here, is only interested in the partition function, then one can 
recast the result of this complicated analysis in terms of the initial variables of BF theory 
subjected to the initial naive primary second class constraints. The only non-triviality that 
remains lies in the B field measure factor. This factor is purely ultralocal, independent of 
the connection and therefore does not affect the analysis or the discretisation of any of the 
constraints to be imposed. The key point about this strategy is that the imposition of the 
constraints is done at the classical level and not at the level of some quantum operator. 

This drastic simplification might not occur if one deals not with the partition function, 
but with the expectation values of some observables, see e.g. [13] for a recent discussion 
of this point. As argued in [23], this ma y s trongly affect the analysis of the operatorial 



version of the constraint as performed in [18|, |23[. However, since in the present work we 



are interested only in the partition function, the strategy is to re-write the reduced phase 
space path integral in terms of the original unconstrained variables, together with a set of 
5-functions imposing all the constraints and gauge-fixings, together with a simple ultralocal 
modification of the path integral measure. 



B. Discretized gravity action and the Barrett-Crane model 

The idea behind the spin foam model approach to quantum gravity is to appropriately 
discretise the gravitational action, and perform all impositions of the constraints, gauge 
fixings, as well as the path integral quantization in this discrete setting. It is very well 
understood how this procedure is realized in the context of topological BF theory. It is then 
the availability of the BF formulation of gravity, which makes gravity, at least superficially, 
look not too distant from BF theory, that makes us hope that one day the gravity case will 
also be understood fully. 

Then our strategy, which is justified by the results of Henneaux et al. [2^ reviwed above, 
is to restrict the measure of integration over the classical configurations that enter into 



the discretisation of the BF path integral. The results of [26| tell us that one just has to 
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impose the primary second-class constraints - the simphcity constraints. As usual, there are 
discretization ambiguities in passing from the continuum to the discrete picture. However, 
as we will see, once we have fixed the discrete set of dynamical variables and have chosen 
on which set of variables to impose the constraints, the implementation of the constraints 
is rather straightforward and almost unambiguous. As we shall see, one ambiguity that 
is left is in the choice of the edge amplitudes and amounts to the choice of the ultralocal 
measure on the B field. Another subtlety lies in the fact that we impose the constraint in a 
linear manner, which involves an introduction of an additional normal vector ni, see below, 
and the full constraint analysis of this case has not yet been shown to be equivalent to the 



one of |26| . In the first archive version of the present work we have given strong evidence 
that the two formulation of the constraints lead to the same measure (we have shown it for 
only half the constraints). For simplicity and coherence we have left this analysis out of 
this printed version, as it is not directly relevant to the main line of argument of this work. 
Another possible worry about our procedure is whether the canonical construction of the 
path integral measure in the continuum commutes with the discretisation. Thus, a canonical 
analysis similar to j26| should definitely be attempted at the discrete level. In spite of these 
shortcomings, the procedure that we propose in this paper seems to be well motivated by 
the continuum results, and, we hope, provides an important new step in the direction of the 
completion of the spin foam quantization programme. 

The above discussion of the ideological foundations of the spin foam method makes us 
ready to dive into more technical intricacies. Thus, when discretizing the theory on a 
simplicial complex one integrates the two-form B^'^ over each face of the triangulation and 
associates to any face / of the triangulation a Lie-algebra quantity see e.g. 0]. One 
also assigns a group element to each tetrahedron and defines the following discrete analog 
of the gravitational action 

S = J2^v{XfGf), (3) 

/ 

where G/ is the holonomy around /. Analogously to the continuous formulation one has 
to impose the "simplicity" constraints on the bivectors Xf, as well as other second-class 
together with first-class constraints on all the discrete variables Xf,Gf. However, in this 
paper we will only concentrate on the simplicity constraints. Some preliminary ideas on 
the other, not considered in this paper second-class constraints, as well as the path integral 
measure are contained in the arxiv version of the present work. However, because of their 
rather inconclusive character we have decided not to include that material in this printed 
version. 

The simplicity constraints are easiest to understand by thinking about a geometric 4- 
simplex in (or Minkowski space M^'^). Given such a simplex one can then associate 
to each face an "area" bivector that we will call Ay that is obtained by taking the wedge 
product of two edge vectors. This way one obtains four bivectors A('^ per each tetrahedron 
t. There are three constraints that this set of bivectors must satisfy: 

1. Simplicity of each face bivectors : 

AyAfu = 0, (4) 

where A^"^ = {l/2)e^'^^^ Akl is the Hodge dual of a bivector. This constraint implies 
that Ay is "simple", i.e., is of the form A^'^ = u^^v^^ and defines a geometrical plane. 
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Cross simplicity : 

Ay A J J = (5) 

for two different faces /, /' in t. These constraints guarantee tliat from tlie four vectors 
defined by two simple bivectors A^'^ only 3 are linearly independent and that therefore 
the planes defined by them span a three-dimensional subspace in M^. 

Closure : 

E4' = 0' (6) 

fct 

This constraint fixes the relative coefficient expressing the linear dependence among 
the bivectors, which insures that the norms of these bivectors can be interpreted as 
of the faces of a geometric simplex. 



For more details on these constraints, see for example [2l[. Importantly, the above con- 
straints are not enough to guarantee that the quadruple of bivectors is one coming from a 
geometric tetrahedron, as there are certain discrete ambiguities. This point is also empha- 



sized in [2l| and plays an important role in the arguments of this paper. At the quantum 
level the bivectors Ay are promoted to the Lorentz group Lie algebra elements. The above 
"simplicity" constraints then have to be imposed on these operators. Barrett and Crane [sl 
realized that the "simplicity" constraints can be imposed at the quantum level by select- 
ing the so-called simple representations of the Lorentz group, intersection constraints are 
imposed by selecting the intertwiner to be the so-called Barrett-Crane intertwiner, and the 
closure constraints can be imposed by selecting only the gauge-invariant intertwiners. From 
the point of view of the discrete action the variables Xy that appear in the action must 

be identified with the Hodge dual Ay of the area bivector A^-^. Below we shall see how this 
procedure can be implemented. 



C. A new geometric criterion for (cross-) simplicity 

In this work we will base our discussion on a somewhat different interpretation of the 
simplicity constraints. This new interpretation will in particular help us to deal with the 
discrete ambiguities that arise. The simplicity constraints, as reviewed above, are quadratic 
in the bivector variables, and this is ultimately the cause of the discrete ambiguities that the 
set of solutions possesses. The new formulation of these constraints makes them linear in 
the bivectors, thus removing the discrete ambiguity. Thus, the new geometric way to look 
at the simplicity and intersection constraints is contained in the following two lemmas: 

Lemma II.l. A bivector X^-^ m (M^^^) is an anti- symmetrized product of two vectors 
if and only if there exists a vector such that X^^nj = 0. 

Proof. If the bivector is simple it defines a two-plane in (M^'^). Taking to be any 
of the vectors orthogonal to this plane proves the assertion in one direction. To prove it in 
the opposite direction we assume that X^"^ is not simple. Then there always exist 4 vectors 
Ui 2, vi 2 so that 

x^-^ = 4V + 4\^ (7) 

and Ui^2, Vi^2 span (M ^'^). Therefore, these vectors can be used as a basis, and the vector 
can be represented as a linear combination of these basis vectors. It is then easy to see 
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that the condition X^'^rij = says that the 4 vectors Ui^2i "^1,2 are hnearly dependent. This 
contradicts our assumption that X^'^ is not simple. □ 

Lemma II. 2. Two simple bivectors X^"^ and X|"^ span a 3-dimensional subspace of 
(M^'^) if and only if there exists a vector such that Xf'^nj = and Xj'^nj = 0. 

Proof. If two bivectors are simple and the 2-planes defined by them span a 3-dimensional 
subspace, then can be chosen to be a vector orthogonal to this subspace. This proves 
the assertion in one direction. The proof in the opposite direction is similar to the proof 
in the above lemma in that one shows that the condition of existence of implies that 
the 4 vectors defined by the 2 two-planes must be linearly dependent, and thus span a 
3-dimensional subspace. □ 

Importantly, the above simplicity and intersection constraints can also be formulated in 
exactly the same way in terms of the Hodge duals X^"^ = {l/2)e^'^^^XKL of the bivectors. 
Thus, we have the following two "dual" lemmas: 

Lemma II. 3. A bivector X^'^ in (M^'^) is is an anti- symmetrized product of two vectors 
if and only if there exists a vector such that X^^nj = 0. 

Proof. A bivector X^"' is simple if and only if its dual bivector X^"^ is simple, so the first 
lemma above also proves this lemma. □ 

Lemma II. 4. Two simple bivectors Xf"' and X^"^ span q 3-dimensional subspace of 
(M^'^) if and only if there exists a vector such that Xl^uj = and X2'^nj = 0. 

Proof. Two bivectors Xfg span a 3-dimensional subspace if and only if their duals Xf 2 do, 
so the second lemma above proves this lemma as well. □ 

An important point is that, in spite of having these two seemingly equivalent ways to 
impose the simplicity and intersection constraints, they are no longer equivalent when one 
starts to deal with more than two bivectors. Thus, considering, as in the previous sub- 
section, the area bivectors Ay obtained by taking the wedge product of the edge vectors 
corresponding to a tetrahedron one can see that Ay satisfy the simplicity and intersection 
constraints in their first version, namely that there exists a vector orthogonal to all Ay , 
but not in the second version. Indeed, let us consider a 3-dimensional tetrahedron built on 
unit vectors along the axes 1, 2, 3. The corresponding area bivectors are: A{2 = e{ A 63 , 
^23 = ^2 A 63 , Al( = e| A e(, and the fourth bivector being — ^{2 — ^23 ~ ^31 • Here 
^0 1 2 3 are unit vectors in the direction of the axes. All 4 bivectors are orthogonal to the 
vector Cq. However, the dual bivectors A{2 = Cq A 63 , ^23 = Cq A e(, A^^( = eg A 63 do not 
satisfy the simplicity-intersection constraints in their first version, as there is no vector such 
that A{2, ^23 ' ^31 would all be orthogonal to. Instead, these bivectors satisfy the simplicity 
constraints in their second version, they are all parallel to cq. 

It is therefore clear that the simplicity and cross-simplicity constraints can be imposed 
on a set of 4 Lie-algebra valued variables Xy in two "dual" ways. Indeed, the constraints 
become the requirement that there exist a vector such that the 4 variables Xy are either 

all orthogonal Xyutj = to the vector or are all parallel Xyutj = to it. This is the 
form of the simplicity constraints that we will use. 

In the first case Xyritj = the variables Xy can be interpreted as the area bivectors 
of a geometric tetrahedron: Xy = Ay , while in the second case Xyutj = the duals Xy 
admit the same interpretation: Xy = Ay . 
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D. Self- and anti-self-dual decomposition 

Consider now one of the variables X^'^ . Consider the following two quantities: 

X^\X) := (niX'J ± v^n/X") . (8) 

Here is at this stage arbitrary unit vector (future timelike in the Lorentzian case) and a is 
1 in the Euclidean signature and —1 in the Lorentzian case. It is easy to see that X^^ni = 0, 
so both vectors X"^ actually take values in the 3-dimensional subspace orthogonal to . It 
is also easy to check that 

X^'{X) =±,/^X^\X), (9) 

where ~ is the Hodge duality operator. Thus, are eigenvectors of this operator. They 
are called the self- and anti-self-dual parts of X^"^ . Their knowledge is equivalent to the 
knowledge of X^'^ (there is an explicit reconstruction formula that we will not need here) . 

As we have seen above, when X^'^ is simple there exists a vector such that X^'^nj = 
or, equivalently, a different vector such that e^'^^^njXxL = 0. The two distinct criteria of 
simplicity (and intersection) can then be formulated in terms of self- and anti-self-dual parts 
of bivectors as: 

I : X^-^nj = X+-^ = X"-^ = n/X", (10) 

n : e^J^^njXKL = ^ X+-^ = -X"^ = v^n/X". (11) 



E. Identification of bivectors with Lie algebra elements 

From now on we restrict our discussion to the Euclidean case. As we have already said, 
the idea of quantization is to first associate to each bivector a Lie algebra element 

X = X^'Jij, (12) 

where J/j are the generators of 5*0(4) 

[Jij, Jkl] = iSixJjL - iSjK-JiL - iSiL-JjK + iSjlJik- (13) 

Let us now pick a fixed four vector uq = (1,0,0,0) and define the corresponding self-dual 
operators: X±^ = (X^ ± X°^), where X* = l/2e'jkX^'' and = (J^ ± Joi). The Lie algebra 
element then decomposes in terms of its self- and anti-self-dual parts 

X = X+ + X_, X± = X^Vf . (14) 

The signs entering the definition of the self-dual operators are such that they form an 
su{2) © su{2) Lie algebra 

[Jf,jf]=226,/J±, [J+,J-]=0. (15) 

At the group level the described splitting of the Lie algebra into the self- and anti-self-dual 
parts leads to the isomorphism S0(4) = SU(2) x SU(2)/Z2, so that each S0(4) group element 
g can be represented as g = {g~^,g~),g'^ G SU(2). This isomorphism can be described more 
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explicitly as follows. A point in ~ can be parametrized by a pseudo-unitary 2x2 
matrix 

Transformations 

SU(2) X SU(2) 3 {g+, g~) : x ^ 9^9^^' (17) 

preserve the determinant det(x) = and are thus orthogonal transformations of M^. 

Of particular interest for us is the "diagonal" subgroup of SU(2) x SU(2) given by elements 
of the form {9,9). Let us also mention that under the adjoint action of 5*0(4) the self- and 
ant i- self- dual elements transform as SU (2) Lie algebra elements: 

gXg-i = 9^X+{g+)-' + g-X^ig-yK (18) 



F. Gravitational and topological sectors 

Now, the key point for us is that, as a comparison between the continuous ([1]) and discrete 
([3]) actions suggests, if the theory is to correspond to the gravitational sector the quantities 
Xy in ([3]) should be identified with the duals of the area bivectors Ay-. 

Xy = le^^KLAf^. (19) 

Correspondingly, the simplicity and intersection constraints should be imposed in their ver- 
sion flTT]) . The other case, namely with Xy satisfying the constraints in their version f[TU|) 
corresponds to the topological sector. In the subsequent sections we will see that there exist 
two corresponding models, one for the gravitational and one for the topological sector. 

Let us rewrite the criteria (fTOj) . (fTTj) in a more convenient form. We will make all state- 
ments about the gravitational sector case. The other case is obtained by a simple duality 
transformation. Thus, as we have just seen, given four face bivectors Xf associated with 
a tetrahedron t and satisfying the simplicity and cross simplicity (intersection) conditions 
there exist a unit vector rit such that nnXy = 0. There always exists an S0(4) transforma- 
tion that brings the vector nf into the fixed reference vector Uq. Let us denote the inverse 
of this transform by {g~^, g~). Thus, no = 1 — >■ n = g^{g~)~^ ■ We have seen that in the case 
when is chosen to point along the "time" axis = (1, 0, 0, 0) the self-and ant i- self- dual 
part of a bivector satisfying the simplicity conditions in their form ffTTl) satisfy Xj' = —Xj 
so that X/ = {Xf, —Xf). The SO (4) rotation maps this bivector into 

{Xf, -Xf) ^ {g-'Xf{g^)~\ -g~Xf{g~)~'). (20) 

This makes it clear that in the case of a general normal vector rit the simplicity-intersection 
constraints on Xf corresponding to the gravitational sector can be formulated as: 

y.f = {Xf,-nyXfnt), (21) 

where Ut G SO (4) should be interpreted as carrying information about the normal vector 
to the tetrahedron containing Xf. Similarly, the constraints appropriate to the topological 
sector take the form: 

^f = {Xf,nyXfm). (22) 
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These are the final criteria we will use in the following. 

In the next section we would like to impose the conditions (12T|) and (l22l) at the quantum 
level. We will achieve this via the technique of the coherent states introduced in the context 
of spin foams by Livine and Speziale [l7| and obtain two distinct models that we will 
interpret as corresponding to the gravitational and topological sectors. The technique of 
coherent states that we shall use will also allow us to see that the model of Barrett and 
Crane does impose the correct gravitational constraints ( 1211) but in a form that is too 
strong: the constraints are imposed for each 4-simplex independently of the others simplices 
of the triangulation. This will be interpreted as a likely source of all the problems the 
Barrett-Crane model is known to suffer from. The model of Engle, Pereira and Rovelli [3] 
on the other hand will be seen to impose the simplicity constraints in a satisfactory manner. 
However, as we will see, this model turns out to give a quantization of the topological sector 
in that it imposes the constraints in their version (l22l) . and is thus unsuitable as a model 
for quantum gravity. A new model that imposes the constraints in their version fl2ip will be 
obtained below. 



III. COHERENT STATES AND THE SPIN FOAM MODEL OF BF THEORY 



The technique of coherent states is crucial for all the constructions that follow. We will 

we will illustrate this technique 



therefore start by a review of this technique. Following [17 
by applying it to the BF theory partition function. 



A. Coherent states 

Let us first describe the coherent states for the group SU(2). It will then be trivial to 
generalize it to SO (4). 

We will start from the following obvious decomposition of the identity operator in the 
representation space of dimension dj = 2j + 1: 

= ^\j^^){j^M^ (23) 

m 

where \j,m),m G [— is the usual orthonormal basis in VK Let us rewrite the same 
operator by inserting in it the Kronecker delta 

Smm' = d, / dg C^. ig)tl,^ig) (24) 

^SU(2) 

written in terms of an integral over the group. Here tl^j{g) is the matrix element of the 
group element g in the representation j computed between the states {j,m\ and \j,j), the 
latter being the highest weight state. More generally: 

tLm'i9) = {j,m\g\j,m'). (25) 

This gives: 

= ^ij] / dgtl^j{g)tl^,Ag) =dj dg\j, g){j, g\, (26) 

^SU(2) 7SU(2) 
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where we have introduced a notation: 

\j,9) = 9\j,j) = J2\^,rn)tl^ig). (27) 

m 

The states \j, g) are known as the coherent states, and the last expression in ( 126|) is the 
decomposition of the identity in terms of these coherent states. 

The decomposition of the identity (!26|) can be further simphfied by noticing that the 
matrix elements tl^j{g) and tl^j{gh) differ only by a phase for any group element h from 
the U(l) subgroup of SU(2). Therefore, the integral in (!26l) can be taken over the coset 
G/H, G = SU(2), H = U(l). Thus, we have: 

lj = dj [ dn\j,n){j,n\. (28) 
Jg/h 

We will suppress the domain of integration G/H in what follows. 

The states \j,n) form (an over-complete) basis in . One of the main motivations for 
using these states comes from the fact that they have a clear geometrical interpretation: Let 
J* be the generators of the SU{2) Lie algebra [J\ J^] = ie^^^Jk- We can easily compute the 
expectation value of these generators in a given coherent state: 

^{j,n)(^i = 'n\T\i, n)ai = j nasn'^ = j n'ai, (29) 

where the quantities CTj are the Pauli matrices and we have denoted by the same letter the 
SU (2) matrix n that parametrized the coherent state and the three-dimensional unit vector 
n that this matrix corresponds to. Thus we see that | j, n) describes a vector in of length 
j and of direction determined by the action of n on a unit reference vector. It is convenient 
to denote this vector by the same letter n. This geometric interpretation will be extremely 
useful for us in what follows. 

An important remark is as follows. Note that we could in principle define a notion of 
"coherent" state starting from any state \j,m) instead of the highest weight vector 
Indeed, we would have the same decomposition of unity as in (l26l) . The geometrical inter- 
pretation would be, however, drastically different. Namely, we could no longer interpret the 
states g\j,m) as representing a vector lying on the sphere of radius j. Instead, these states 
correspond to a circle on this sphere of radius y'j^ — m^. Moreover, we can compute the 
dispersion of the quadratic operator with the result being 

AJ"^ = J +f -m^. (30) 

Thus, we see that only the highest and lowest states m = ±j lead to the coherent states 
that minimize the uncertainty relation. 

It will be of crucial importance for us that there is not one but two coherent states 
g\j,j) = \j,g) and g\j, —j) = \j,g). In order to relate them let us remark that the group 
element 

1 
-1 



(31) 



is such that 

ge = eg, (32) 
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where g is the complex conjugate matrix. The matrix element of e in any representation is 
given by 

= (-l)^""5n,-,n. (33) 

Thus, we have (j, — j|e|j, j) = (—1)^'' and therefore the lowest weight vector can be obtained 
by the action of e on We have: 

\j,-j) = (-lf^e\j,j). (34) 

More generally we can define the conjugate states 

\j-^ = e\j,m) = (-l)^+™|j,-m). (35) 

These are called conjugate since their matrix elements are the complex conjugates of the 
usual matrix elements 

JlMgM = {j,m\e~'ge\j,n) = = (-l)'^+"+"t-„^-n(^?) (36) 

We thus obtain two families of coherent states \j, g) and 

^ {-lf^g\3, -j) = Y.lJ^'^)- (37) 

m 

We are eventually interested in coherent states for SO (4). What this discussion shows is 
that one has four possible coherent states which are given by acting with an SO (4) group 
element on either of the following states: 

|J,-J) ® ® (38) 

Note that these four states can be obtained from one another by the action of an 5*0(4) 
group element. However, if one considers only the action of the diagonal SU(2) subgroup 
of elements of the form {g,g) then there are two inequivalent states that cannot be related 
by such a transformation. As we will see below the state \j, g) ® |j, g) will play the key role 
in the constructions of the topological sector whereas the state | j, g) ® |j, g) will be the one 
relevant to the gravitational sector. 

B. BF theory 

In the approach to spin foam models that derives them from an action principle, the 
starting point is the fact that gravity can be written as a constrained BF theory. Thus, 
the strategy devised in [sl is to first construct the partition function of the BF theory on a 
given triangulation and then impose the simplicity constraints at the level of the partition 
function. Note that this is rather different from the strategy followed by Barret and Crane 
who imposed the simplicity constraints at the level of each 4-simplex, with the starting 
point being the idea of quantization of a 4-simplex. In simplicial version of BF theory, the 
partition function is obtained by imposing the flatness condition on a discrete connection 
— using a product of J-functions imposing the constraints that the holonomies around all 
dual faces are trivial. These constraint are obviously too strong for a gravitational theory 
and are to be relaxed. 
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To write these constraints more explicitly, we introduce the following convenient nota- 
tions. Let us consider two neighboring 4-simplices a, a' of a simplicial decomposition of a 
4- manifold. These two simplices intersect along a tetrahedron t = {a, a'). One chooses "cen- 
ters" of these 4-simplices (these centers become vertices of the dual complex), and connects 
them with an interval - the edge aa' of the dual complex. This dual edges are in one-to-one 
correspondence with a tetrahedron t of the original simplicial decomposition. Consider now 
the holonomy gt = Qaa' of the gravitational spin connection along the dual edge aa' . Assum- 
ing that the bundle is trivial (or choosing a local trivialization), this holonomy is an element 
of the group S0(4) in the Euclidean case and SL(2,C) in the Lorentzian setting. Taking a 
product of these group elements for all the dual edges that form the boundary of the dual 
face (one has to introduce an orientation of the dual face in an arbitrary fashion) one obtains 
the holonomy around the dual face. The flatness condition is imposed by requiring all these 
dual face holonomies to be trivial. The partition function of the theory is then given by: 

n^^^riMn^*)' ^^^^ 

t f \tDf J 

The integral here is over group elements associated with each tetrahedron. The (5-functions 
impose the constraint that the holonomy around each dual face is trivial. 

One then uses the Plancherel decomposition to rewrite the 5-functions into a sum over 
irreducible representations of the group in question. In the case of SO (4) we use the notation 
g = ((?+, g-) for an S0(4) group element, j = (j+, j_) for an S0(4) irreducible representation, 
dj = dj_^_dj_, dj = 2j + 1 for the corresponding dimension and trj^(g) = tTj+{g^)tTjj{g~) for 

the characters. It proves instrumental to split each group element gt in (l39l) into a product 
of two group elements. Given that the group elements gt have the interpretation of the 
holonomy along the dual edges, they can be represented as a product of the holonomy g„t 
from the center of 4-simplex a to the center of tetrahedron t with the holonomy from t to 
a'. Thus, gt = gatQtu'- The introduced group elements can be integrated over separately 
without changing the partition function (we assume that the group integration measure is 
normalized). We will also use the convention that gta = g~t- The Plancherel decomposition, 
together with the introduced trick of doubling the group elements to be integrated over, 
reduce the BF partition function to the following expression: 



^ = E n / n dg*<^ n n g^^^g*'^' i • (^o) 



The key step proposed in [17| is to insert into this expression a multiple decomposition of 
the identity written in terms of the coherent states. Thus, for each tetrahedron t and each 
dual face we insert into (140|) the following decomposition of unity: 

Ij = dn|j,n)(j,n|, |j, n) = n_) ® n_). (41) 

This gives 

if f (t.<^) {*,/) i^J) 
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For a triangulation of a 4- dimensional manifold there are exactly 4 dual faces that share each 
particular dual edge, and thus each group element gta enters into 4 characters. Integrating 
over these holonomies one produces a product of 15j symbols - one per 4-simplex - and gets 
what is known as the spin foam model of BF theory j28|. This derivation is usually done 
without the insertion of a complete set of states for each tetrahedron, which only complicates 
the matters in this simple case. However, below we will see that the expression (H2l) that we 
wrote is a very useful starting point for the non-topological models. 

Another useful interpretation that we can give the expression P2|) is to think about 
the triangulated manifold amplitude that it gives as the product of "vertex" amplitudes 
corresponding to 4-simplices a (or vertices of the dual triangulation) and "edge" amplitudes 
corresponding to tetrahedra t (or edges of the dual triangulation). The vertex amplitude is 
a function of the representations {j/} labelling the faces of a as well as 5 x 4 group elements 
Mtf labelling the faces / of tetrahedra t. Note that the same face / is shared by exactly two 
tetrahedra, which we denote by t(/) and t'{f). The vertex amplitude is given by: 

^BF[{j/},{nt/}] = n^j/ / n^stn^j/''^*(/)/i^g*(/))"'st'(/)ij/,nt'(/)/)- (43) 

f t f 

The edge amplitude that corresponds to P2l) is simply a (5-function that sets the variable 
^t{f)f °f '^^^ 4-simplex a to be equal to a similar variable ^t'{f)f corresponding to the same 
face / but in another 4-simplex. Thus, the edge amplitude is: 

^Bpiih}^ {<f}] = '^«(/)/' <''(/)/) or equivalently = dj^, (j/, <''(/)/)• (44) 

The full partition function is then obtained by taking a product of vertex and edge ampli- 
tudes, summing over jj and integrating over n^j. 

IV. IMPOSING THE SIMPLICITY CONSTRAINTS 

In the previous section we have seen that it is possible to associate to each state n^j) 
a bivector 

^(j/,n,,) = (j/> nt/l J"|j/, ntf), X(j^,n,^) = „+ )' ^a7,n- )) ' ('^^) 

where J^"^ denotes a basis of SO (4) Lie algebra elements. According to the results of section 
ini demanding that these bivectors satisfy the geometrical simplicity and cross simplicity 
conditions i), ii) is equivalent to the requirement that there exists an SU(2) group element 
Ut (interpreted as characterizing the normal to the 3-dimensional plane spanned by the 
tetrahedron t) such that 

^ih;n.f) = iXf,-UtXfu;'). (46) 

This condition is the one relevant for the gravitational sector. The topological sector corre- 
sponds to a similar condition, but without the minus sign in the anti-self-dual part. 

Using the fact that Xq-^^^^) = —X{jf,n) and the fact that X(^j^^nh^) = ^(j/,n) if = 
exp(i0(j3) is a U(l) element, we can rewrite the simplicity conditions for the gravitational 
sector in terms of the spins jf and the SU{2) elements rij as: 

Jf=jJ, and {n^j,n;'^) = {ntfh^^f,utntfh^\), (47) 
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where (j)tf G [0, 27r]. The reason why we have allowed for the presence of an additional U(l) 
comes because the group elements parametrizing the normals are defined modulo the 
right U(l). Note that the bivectors are initially invariant under U{1) x U{1). However, the 
diagonal U{1) invariance is reabsorbed into the definition of n/ and the anti-diagonal one 
is taken care by (ptf- We will see below that these angles can play an interesting role in the 
construction. 

Similarly, the simplicity constraints for the topological sector are equivalent to the fol- 
lowing two conditions: 

Jf=jJ=jf, and {nlj,n~f) = {ntfh^^^,utntfh^^^), (48) 

with the only difference being in the absence of the e group element in the anti-self-dual 
copy. 



A. The Barret- Crane model 

The logic now is to impose the constraints fHTl) and psj) in the expression P2|) for the 
partition function of the BF theory, with the motivation being that a similar imposition 
at the level of the continuous theory is known to produce gravity. It is clear that when 
doing this, one will be imposing the constraints at each tetrahedron. Importantly, as the 
edge amplitude 0441) requires, the constraints are not imposed at each 4-simplex individually 
since a tetrahedron is shared by two 4-simplices. However, let us see what happens if we 
ignore this requirement and do impose the constraints at the level of each 4-simplex. Thus, 
what we would like to do is to impose the gravitational constraint (147!) on the group elements 
entering the definition of (l43l) and then integrate over ntf, (ptfiUt- The integral over ut 
can be reabsorbed in a redefinition of and the integral over (f)tf also drops out because of 
the simplicity constraints = j~ ■ It is interesting to note that even if one does not impose 
the constraint = j~ it will be automatically implemented by the integral over the U(l) 
group element (ptf- This gives another reason to introduce this group element - to give a 
simple tool to impose the constraints = j~ . After all this operations are done, one is left 
with the vertex amplitude that only depends on 10 spins jf and is given by: 

= [U^gtdgtU'^n.fX (49) 

/ 



where we have used the fact that \j,ge) = \j,g). The integrals over utj can be easily taken 
as follows: 

m,m' 

= ^ E \j^^)^\j^^'Km'i9^i9-r')- (51) 



m,m' 



Using this result and some obvious changes of variables of integration, it is not hard to reduce 
l9|) to the following, well-known in the literature, expression of the vertex amplitude of the 



17 



Barrett-Crane model: 

AbcUjA] = [U'^at n^^^/((^*(/))"'^?*'(/))> (52) 
t f 

where {jf} is the collection of 10 irreducible representations that label the faces of a, 
t{f),t'{f) are the two tetrahedra of a that share the face /, and tTj{g) is the character 
in j-th representation of an SU(2) element g. We thus see that the obviously wrong way of 
imposing the constraints that ignores the fact that the face normals ritf in two neighbouring 
4-simplices should be the same produces the Barrett-Crane amplitude. The above derivation 
both pinpoints what is wrong about the Barrett-Crane model, as well as suggests a way to 
remedy the problem. 

One very important feature of the Barrett-Crane amplitude, which is also obvious from 
our analysis, is that it clearly corresponds to the gravitational and not the topological sector. 
This fact is in agreement with the analysis of Baez and Barrett [2l| where it was shown, 
using very different arguments, that the Barrett-Crane way of implementing the constraints 
cannot lead to a quantisation of the topological sector. 

It is illuminating to see the impossibility of having the topological sector "Barrett-Crane" 
model via our technology. Indeed, in an analog of the Barrett-Crane model for the topological 
sector, the amplitude should be obtained by implementing the topological constraints 



and therefore using the coherent states |j, ritf) ® |j, ritf) instead of \j, ritj) ® \j, ntf). However 
the integral 

m,m' 

= Y\j,rn) 0\j,m') j dntl^j{g^n)tl^,j{eg^ne~^) (54) 

m,m' 

= ^ E l-^'"^) ®\j^^')tLm'i9^i^9-r%A^) = (55) 

m,m,' 

is vanishing since tjj{e) = 0. Thus, one clearly sees that there is no way to impose the con- 
straints of the topological sector a la Barrett and Crane on each vertex amplitude separately. 

We have not discussed the closure constraints so far. However, these constraints are 
automatically imposed by the integrals over the S0(4) holonomy group elements g^^f Indeed, 
these integrals "average" the sum of all four bivectors corresponding to a tetrahedron over 
all directions. Thus, only the configurations that add up to zero survive the integration. 

The above derivation has made it clear that, despite the important feature of the Barrett- 
Crane model as corresponding to the "right" gravitational sector of the theory, it is, however, 
plagued by a deficiency. Indeed, as we have seen, in order to construct the Barrett-Crane 
model one has to integrate separately over the n^j variables for each 4-simplex. But since 
the variables n^j describe the bivectors corresponding to the faces of a tetrahedron t, the 
Barrett-Crane model describes "geometry" in which a tetrahedron t viewed from one simplex 
has a completely different "shape" as when viewed from a different 4-simplex. This is clearly 
not a geometrically adequate way of implementing the constraints. A related comment is 
that the Barrett-Crane vertex amplitude does not depend on the normals Utf, which means 
that some of the geometric information about the geometry of the triangulated manifold is 



lost. It is possible that the fact [Ij], [16[, [15| that this vertex amplitude is dominated by 



degenerate tetrahedra is precisely due to this loss of geometrical information. 
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B. The topological sector and the model of [1 

From the analysis performed in the previous subsection it became clear that the main 
deficiency of the Barrett-Crane model is easy to resolve. The idea for how to do it is a simple 
extension of the set of ideas developed in jsj. Namely, the continuum partition function of 
BF theory is expressed as an integral over the connection field A and the 2-form field B. In 
the previous section we saw that at the discrete level the BF theory partition function can 
be expressed as a double integral over the connection (holonomies) and face variables X^j, 
see ( 142|) . where X^j are described by the group elements n^j. Similarly to the continuum 
theory, one should impose the simplicity constraint at the level of the measure of integration 
over X(/. Importantly, since each tetrahedron t is shared by two 4-simplices a, a' the same 
variable X^j does appear in two vertex amplitudes. Imposing the constraints at the level of 
the integration measure over X^j means that while we ask these bivectors to be simple we 
also impose the condition that X^j = X^j. Geometrically this corresponds to the condition 
that the geometry of the tetrahedron t viewed from the point of view of 4-simplex a agrees 
with the geometry viewed from a'. 

We start by applying the above considerations to the topological sector of the theory. The 
constraints we have to impose are: (i) instead of summing over all SO (4) representations j/ 
we sum only over the simple representations j/ = (j/, j/) and (ii) instead of integrating over 
all group elements ntj we integrate over the normals of the form 

After some simple changes of variables the integration the integrals over Ut and h^^^ can be 
dropped and we are left with the following amplitude: 



n (Jf^'^tflidat)"^ 9at'\jf^i^t'f){jf,ntf\{g-t)-^ 9at'\jf^i^t'f) (56) 

where all the integrals are over SU(2) elements. 

As in the case of the Barrett-Crane model it might seem that only two out of three 
constraints - namely, simplicity and cross-simplicity, but not the closure were imposed. 
However, as before, the integration over variables g^t does impose the closure constraints 
since 

J2xi-Uf I <lg\3f,ghf)\=Q. (57) 

It is clear from the expression above that the essential difference between this model and 
the model for BF theory is in the form of the intertwiner which is inserted for each pair ft. 
Indeed, instead of inserting the decomposition of the identity of the BF model: 



lj+ = dj-dj+ j dn dri^ \i ,n )®|j'^,n+)(j ,n (58) 
the model we are discussing is obtained by inserting the simple intertwiner 

Tj =d) j dn I j, n) ® | j, n) {j, (g) (j, n| . (59) 
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One can compute this intertwiner by recalling the standard tensor product isomorphism 
Vj ^Vj = (sfto^k^ which is realised in terms of the group invariant maps. Thus, there exist 
invariant maps 



C^j k . Yj+ ^ yj _ yk^ fjj+j^k . yk _^ yj^ ^ yj- _ (gg) 

These maps are unique up to normalisation, we chose the normalisation such that 

j + r 

J2 dkC^^^-'^C^^'''' = lj+^lj~, C^^rk(jj+j~k' = J^i^_ (61) 

k=\j+-j-\ *^ 

The matrix elements of these intertwiners are the Clebsch-Gordan coefficients (3j symbols) 
and their complex conjugates 



(62) 

It is now not hard to show that 

® |j,n) = ^C^(|2j»), (63) 

where we have denoted = C^^^K The fact fl63|) is of importance, so it is worth seeing how 
this comes about in detail. Let us recall the definition (|27|) of the coherent states in terms 
of the matrix elements. Using this we can write: 

|j, n) ® |j, n) = ^ |j, m") ® |j, m'^)t^^_.{n)t^^+.{n). (64) 

m~ ,m.+ 

We now use the formula fll54p for the product of two matrix elements involving the Clebsch- 
Gordan coefficients to re- write: 

2j 



2+ ,m,m' 1=0 



However, it is easy to see that the last Clebsch-Gordan coefficient is only non-zero for 
m' = j + j = 2j, and that the only representation for which this is possible is k = 2j. The 
Clebsch-Gordan coefficient of relevance is computed in the appendix, with the result being: 



^S' = (66) 



This implies that 



|j, n) ® |j, n) = Yl l-^' "^") ® l-^' ^^)Ci~m+Aj{n) (67) 

= v^^C^(|2j,m))t^2,.(n) = v^C^(|2j,n)), 

m 

which proves (1631) . We can now easily perform the integral 

Tj = d^^daj J dna\2 j, n){2j,n\ = d^j & o ^ . (68) 



20 



Note that this intertwiner is proportional to a projector TjTj = 5j j(d^/d2j)Tj. 

This intertwiner is (up to normalisation) the one introduced by Engle, Pereira Rovelli in 
(l8| . This shows that the model we obtain after integrations over the group elements g^^ 
is (up to a different edge normalisation) the model of [l8|. However, it is clear from our 
analysis that this model gives a quantisation of the topological, not the gravitational sector. 
Thus, despite providing an important step in the right direction, the model of l3] is unlikely 
to give an interesting model of 4-dimensional quantum gravity. 

For completeness of the exposition, let us show explicitly how one recovers the model of 



18( 1 as it is presented in this reference. First, let us introduce the following intertwiner: 

y,(ji, . . . , J4) = ^ C^u-^^lz, m) ® C^^=^'\i~^ -.Vj,®---® Vj, C. (69) 

This intertwiner appears in the result of evaluation of the following group integral 

I dg t^^ig) ® ■ ■ ■ ® f'{g) = 5^d,r,(ji, ■ ■ ■,u)Yiij,, ■ ■ (70) 

The problem of computation of the partition function (156|1 of the topological spin foam 
model reduces to that of evaluation of the integral 

I = J dg^dg- ®ti {{t''i9^) ® t^'ig-)) C^\^J^, n,)) . (71) 
This can be done by a repeated use of (1701) . One finds: 

1= Yl d,+ d,- IV (ji) ® W{ji)ftu- Ui) {dkYk{2ji) ®ti |2ji, Ui)) , (72) 

where 

f^+^~ (j.) = {Yi^ Ui) ® Y,- (j,) I ®tl CJ^lnm) (73) 

are precisely the coefficients introduced in [l8|. The vertex amplitude is then given by: 



The dimension factors we obtain are not all included in T^. Apart from that, the amplitude 
is exactly the one proposed in this reference. The final model is obtained by summing over 
ke and jf with the measures d^ and ^ • d^ and by associating to each strand of the dual 



edge (equivalently each pair tf) a factor of d^. This is the model described in [l8| with 
a slightly different choice of the edge amplitudes. Note that the precise form of the edge 
amplitude depends on how the simplicity constraints are imposed. 

The model of [l8| thus solves the main problem of the Barrett-Crane model: it does 
solder the geometry of one tetrahedron to the geometry of the same tetrahedron viewed from 
another 4-simplex. That in itself is a great success. The use of the coherent states renders 
the construction leading to this model rather transparent in that it recovers it from simple 
and well-motivated operations on the partition function of the BF theory. However, as we 
have already emphasized above, the model described gives a quantisation of the topological 
sector of the constrained BF theory. 
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C. The gravitational sector: a new model 

The two models described above, namely that of Barrett-Crane and that of Engle et al 
each has its own merit, but, as was clear from our analysis, both fail to provide a satisfactory 
model for 4d quantum gravity. Our analysis however also suggests a way to merge together 
the advantages of the Barret-Crane model (as corresponding to the gravitational sector) 
with the advantages of Engle et al. model (soldering the neighboring geometry). This leads 
to a new spin foam model free of the described drawbacks. 

The idea is essentially the same as the one already implemented in the case of the topo- 
logical sector. Thus, we impose the following two constraints on the partition function of 
the BF theory: (i) only the simple representations will be included in the state sum and 
(ii) instead of integrating over all SO (4) group elements n^j we integrate only over the one 
having the form rijj = {ntfh^^^.,Utntfeh'^^^). As we have seen above the effect of the group 

element e is to select the coherent states |j, n) ® |j, n) instead of \j^n) ® \j,n). After some 
simple changes of variables of integration the integrals over ut and h^^^ are again dropped 
and we are left with the following partition function: 



'vGrav 



5:ndi/n4d«'//nd9id,„ 

is f (tJ) M 

Y[{jf,ntf\{g+y^ gtt'\jf,nt'f){jf,ntf\{g^t)-^ g^t,\jf,nt'f), (75) 

where all the integrals are over SU{2) elements. 

We therefore see that imposing the simplicity and cross-simplicity constraints as appro- 
priate for the gravitational sector means that the following intertwiner is inserted into each 
"strand" (corresponding to a face /) of each dual edge (corresponding to a tetrahedron t) 



= j n) ® I j, n) {j, n\(g) {j,n\. (76) 

Using arguments similar to those in the previous section one can show that 

\j,n) ® \j,n) = Y,^,^Cia^^ {t\n)%Q)) , (77) 

A:=0 

where the Clebsch-Gordan coefficient 

^k-\y,^jo) - (2j-A;)!(2j + fc + l)! ^ 

is computed in the appendix and |fc,0) denotes the SU(2) invariant (spherical) state in the 
representation V^. 

To prove ( 177j) we consider the tensor product 

|J,J)®|J,-J) = Y.'^kCJrkoC"W3,j)®\j,-3)) (79) 

k 

= J2<^kCJi^\k,m)q%, = Y,dkCJrk\k,o)qi%, (so) 

k,m k 
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where in the first hne we have used the decomposition of the identity in terms of the Clebsch- 
Gordan coefficients, and in the second line the fact that the relevant Clebsch-Gordan coeffi- 
cient is zero unless m = 0. This establishes the result (177|) for the case n = Id. The general 
case follows by considering the action on this state of the matrix P{g) ® t^{g) and using the 
intertwining property of the Clebsch-Gordan coefficients. Using (1771) one obtains: 

k=0 

We note that for k = this intertwiner is the product of the usual Barrett-Crane intertwiners. 
Moreover, as the coefficients Cl decrease with k increasing 

the Barrett-Crane intertwiner is in a certain sense the dominant term in the sum (IHT]) . 

To finish with the presentation of the model we describe it explicitly in terms of the 15j 
symbols. Similarly to what we have done in the previous subsection, we integrate over the 
SO (4) holonomy group elements. The result is most conveniently presented by introducing 
the set of coefficients: 



k) ^ ® Y,-U^)\ ®ti Cnnk.\Y,{k,)). (83) 

In terms of these the vertex amplitude is given by: 

The full model is obtained by summing over kt, jf and ktf with the measures 

En*. En^-.. En^ct- m 



Let us note that, similar to what happens in the case of the model [18|, the boundary spin 
networks of the model described are 4-valent. The S0(3) spins ktf label the edges while It 
label the vertices of the boundary SO (3) spin networks. 

V. THE IMMIRZI PARAMETER 

An important quantity that figures prominently in the loop approach to quantum gravity 
is the Immirzi parameter, customarily denoted by 7. The importance of this parameter 
is illustrated by the fact that it enters in the combination HG'-f into the expressions for 
the spectra of geometrical operators, see e.g. j29|. At the classical level this parameter is 
present in the classical action, and, even though does not modify the equations of motion 
of the theory, it does modify its symplectic structure. This parameter arises in our context 
as being present in the relation between the Lie algebra valued two-form field B^'^ of the 
BF formulation of gravity and the frame field one-forms e^. Thus, we have, in the case of 



non-trivial 7, see [30 



i?" = i^e^'KL + ^5^X^ ) A e^. (86) 
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One of the puzzles associated with the Barret-Crane model has always been the fact that it 
was not possible to construct a spin foam model incorporating the Immirzi parameter, see 



3l| for a discussion and some early developments on this problem. 

We would now like to describe how the approach we have developed in the present paper 
sheds new light on this problem. As we have already seen above, at the discrete (simplicial) 
level an important role is played by the relation between the Lie-algebra-valued variables 
Xy and the " area" bivectors Ay , see (fT9l) . Introduction of the Immirzi parameter changes 
this relation to: 

xy = Ay + ^Ay ^ a^^ = j^ix^^ - -fX^y, (s?) 

where we have assumed that a = 1 for simplicity (Riemannian signature). The original 
gravitational case f|T9l) can be recovered from this general relation by sending j —>■ oo. 

We have already seen that the geometrically correct way to impose the simplicity and 
cross-simplicity constraints is to demand that there is a vector nj such that all the area 
bivectors Ay are orthogonal to . In the case of a non-zero Immirzi parameter this leads 
to the requirement: 

njA^'^ = ^ njXy - -fnjXy = 0. (88) 

When written in terms of the self-dual decomposition this is rephrased as the requirement 
that there exist X in 5u(2) and u G SU(2) such that 

^f=((^ + -) Xf. uXfuA . (89) 



It is now clear that there are two distinct cases to consider: 7 > 1 and 7 < 1 (only the 
first of these cases was discussed in the original arxiv version of this paper). For 7 > 1, 
the quantum version of the simplicity constraints amount to the following restriction on the 
irreducible representations of SO (4): 

7 > 1 3^ = ^ J, (90) 

7 

where j^, j are half integers. It is clear that this condition can be satisfied only if 7 is 
rational, and given by 

7 = > 1- (91) 

3^ -r 



Moreover, in the case 7 > 1 the coherent states to be used when imposing the constraints 
3 the ones in their " 
nus sign in front of 
For 7 < 1 we have: 



are the ones in their "anti-parallel" version |j, g) ® |j, (?), as is clear from the presence of the 
minus sign in front of the second entry in 



7 < 1 j± = 13. (92) 

7 



The corresponding rational value of 7 is: 



7^^<l. (93) 



It is clear that in this case the coherent states to be used are the "parallel" ones |j, g) ® |j, g). 
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From (19T]) . (p3|) we see that the Immirzi parameter is quantized. The case 7 = oo is the 
one previously studied where = j~. Now, for either of the cases (19T1) or (p3!) let us write 
the irreducible representations of SO (4) compatible with the simplicity condition as: 



= (94) 



where j G Z/2 is a half integer and 7^ are as follows from ( !9Ql) or (l92l) . 

The above considerations allow us to write a spin foam model that corresponds to a 
non-zero 7. There are two different models, depending on whether 7 is greater or smaller 
than unity. As before, one needs to specify only the intertwiner associated with each pair 
tf. For the 7 < 1 case the appropriate intertwiner is given by: 



7<1 = dj^+dj^- / dn |j7+,ra) ® |j7 , n) (j7+, ra| ® (j7 ,n\ (95) 



= dj^+dj^- CJ-7+J7-i(7++7-)C^7+i7 Jh++J )_ (96) 
In the 7 > 1 case the intertwiner is: 



7 > 1 G] = dj^+dj^~ J dn |7+j,n) (g) I7 j,n)(7+j,n| ® (7 j,n| (97) 

i(7++7~) 



= d,-^+d,-^- y dfc C^."'^'."'~J, + . ' C'^+J-,7-J,fcC'^-7+i7-fc_ (98) 

J7 J7 / ^ J7+J7 j(7+-7 ) ^ ^ 

k=j{-y+--y-) 

The relevant Clebsch-Gordan coefficients are computed in the appendix. 

One sees that both T'^ and do not contain the trivial representation in their de- 
composition. This gives one explanation of why it was impossible to find them using the 
Barret-Crane way of imposing the constraints. 

One specially interesting case is that of 7 = 1, or, equivalently (7"'' = 1,7~ = 0). In 
this case one might naively expect to get a self-dual model. However, as it can be seen 
from (IHTI) . this limit is singular and the geometrical interpretation of Xf is lost. Thus, the 
self-dual formulation of gravity should not be expected to arise from this limit, at least not 
in a naive way. We note that, interestingly, for this value of 7 the intertwiners of both the 
gravitational and topological sectors simplify and are equal to the identity operator 

a] =T} = 1,. (99) 



VI. A NEW MODEL: LORENTZIAN SIGNATURE 

An added benefit of our construction is that it is possible to generalize the whole discussion 
to the Lorentzian case. 

In case of the Lorentz group we have an isomorphism S0(1, 3) ~ SL(2, C). This isomor- 
phism can be realized explicitly by considering the action of SL(2,C) on Hermitian 2x2 
matrices: SL(2,C) 3 g : x ^ gxg\ where g'^ = (g*)^ is the Hermitian conjugation. This 
transformation preserves the determinant of x, and is thus an orthogonal transformation. 
As in the Euclidean case, there is the "canonical" SU(2) subgroup, which is just the natural 
SU(2) subgroup inside SL(2,C) and which fixes the vector in M^'^ that is represented by 
the identity matrix x. The Lorentzian Barrett-Crane model is obtained by integrating over 
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precisely the described diagonal SU(2) subgroup in each strand forming the dual edge. As 
we will see, a more illuminating way to derive the Lorentzian Barrett-Crane model is via 
the coherent states. 

The key question for us is which coherent states should one work with? Coherent states 
for the Lorentz group have been considered in the literature 35, 3^, but the context of 



these references is not quite what we need. One possible notion is to start with an SU(2) 
invariant state |p, 0) and act on it with a group element to define |p, (7). Because the state 
|p, 0) is SU(2)-invariant, this state is labelled by a point in H^. However, this is analogous 
to considering the coherent states arising from | j, 0) instead of \j,j) in the case of SU(2), 
which is not the correct prescription, as the uncertainty is not minimized in this case. As 
we shall see below, the correct coherent states turn out to be labelled by a point on the 
(projective) null cone. The following subsections develop the corresponding techniques in 
some detail. Once these states are identified the construction of the Lorentzian gravitational 
model follows the same lines as in the Euclidean case. 



A. SL(2, C) simple representations. 

Let us first review how the irreducible unitary representations of the Lorentz group are 
constructed. According to the general representations theory of non-compact groups devel- 



oped in e.g. [3J], irreducible representations are easiest to describe in the space of functions 
on the so-called space of horospheres in the group, which is the homogeneous group manifold 
N\G, where N is the largest nilpotent subgroup of G. For G = SL(2, C) the nilpotent sub- 
group can be chosen to be the group A^_|. of upper diagonal matrices with 1 on the diagonal. 
Let us parametrize the space B := A^+\SL(2, C) by two complex numbers Zi, Z2 G C*, where 
C* is the extended complex plane. The parametrization is as follows: 

b{z,w) = ( ° y (100) 



(101) 



Zl Z2 ^ 

The right action of the group G on B translates then to the so-called affine action 



^'(^1, ^2) ^ J ^ ~ b{(^zi + CZ2, bzi + dz2) 



where the equivalence is that of the left multiplication by an group element. The 
simplicity of the above action of G on i3 allows for an immediate description of the irreducible 
representations: they are realized in the spaces of homogeneous functions on B of some 
degree: 

F(A^i, A22) = A^A^'F(^i, Z2), A, (T, a' e C. (102) 

Another, often more convenient realization of these representations is in the space of func- 
tions of one complex variable, which is possible in view of the above homogeneity property. 
Thus, we introduce a new complex variable z = Zi/ Z2 on which the group acts by fractional 
linear transformations. We then have: 

F(zi, Z2) = z^zY'fiz), f{z) = F{z, 1). (103) 

The transformed function is, on the other hand 

Tg o F{zi, Z2) = {hzi + dz2Y{hzi + dz2f F (^ 11'^ ^2l ' ^) ^2 ^"2"'Tg o Fiz, 1), (104) 
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which defines the action in the space of functions of one variable: 



t;'- o f{z) = {hz + dYihz + d) f[ ^—^ ) . (105) 



az 



A detailed analysis of S^l shows that the unitary irreducible representations are parametrized 



by a pair {s, p), s E Z, p E M+ and correspond to: 

a = —s — ip — l,a' = s — ip — 1. (106) 
The group action is therefore: 



2s 

I n ^ r \ 

(107) 



and the Hermitian product simply 

{f\g) = - I d'zf{z)g{z). (108) 



TT 

The series of unitary irreducible representations (0,p) and (s,0) are often referred to as 
simple. Loosely speaking, the two series can be referred to as simple "continuous" and 
"discrete" correspondingly. Note, however, that there is no discrete series of representations 
of SL(2,C) of the type familiar from the context of the real group such as SL(2,M). The 
simple representations play a special role for they admit a nice geometric interpretation in 
terms of simple bivectors. Namely, let us compute the two Casimirs. Denoting by J^^ the 
Hermitian generators of SO (3, 1) these are given by 

^ J^^J^^ = .2 - p2 - 1, \j^^JAB = 2sp. (109) 

Therefore, we see that s = or p = solves the simplicity constraints and that s = is the 
sector where J is timelike hence J spacelike. 

There exists another convenient realization of representations (s, p) in terms of functions 
(for the case of s = 0) or sections of some bundle (for s 7^ 0) over the null cone or over the 
hyperboloid H^. The case of representations (0,p) is the simplest one and we will restrict 
our attention to it. 

From now on we will use the Dirac bra-ket notation and denote f{z) = {f\p, z). As we 
will see the states \p, z) will be our coherent states in the Lorentzian case. In order to study 
their properties as well as their relation with the Euclidean coherent states it is convenient 
to employ two different realisations of the representations (0,p). 

The first realization is obtained from the fact that SO (3, 1) acts by its geometrical action 
on the space of functions on the null cone. The simple representation (0,p) is obtained by 
considering the functions which are homogeneous of degree minus 1 + zp, that is 0(A^) = 
The intertwining operator between this realisation and the realization fll07p 
is such that the value of this function on the lightcone along its intersection with the plane 
Xo + x^ = 1 exactly matches f{z) (this intersection is the complex plane C*, so this condition 
makes sense). Namely given a complex number we can construct a null vector represented 
as zero determinant Hermitian matrix 
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and the identification is simply (j)f{^{z)) = f{z). In terms of tlie Dirac notation tliis means 
that we introduce states \p,C,) which are such that 



IP' 



z) = \p,a^)), \p,\o = \\r^'''''^\p,o- (111) 



The intersection of the null cone with the the plane Xq = 1 is a two sphere. We can 
therefore parametrise the null vector up to a scale by a point n{z) on the 2-sphere. For 
instance ^{z) = (1 + \z\'^){l,n{z)) where n{z) = n{z)a3n{z)~^ and n{z) is represented by a 
unitary matrix 

^(^) = rr^2 ( --Z I ^ ■ ^112) 



\z\ 



We therefore have a realization in terms of function on or states |p, n), the relation with 
other realization is given in term of the coherent states by 

\p,n{z)) = {l + \z\'Y+^^\p,z). (113) 

The second realization is in terms of functions on H^. The representation (0,p) can be 
obtained from functions on which satisfy the wave equation 

□H30(X) = -(p2 + l)</.(X). (114) 

The scalar product is normalised to be 

(f\9) = J^^-^J{X)9{X). (115) 

The intertwining operator between this and the previous realization is given as follows. 
One can use function on the light cone as the boundary data for such solutions of the wave 
equation f lll4p . It is important to remark that a solution of (11141) can be obtained from 
a solution of the equation □ki,3 0(X) = provided that (f){X) is extended inside the light 
cone as an homogeneous function of degree —{1 + ip) or of degree —{1 — ip). The fact that 
there is these two possibilities is an important subtle point. It is related to the fact that 
the representations (0,p) and (0, — p) are equivalent. More explicitly, the relation between 
the two realizations is as follows. A point x G i^s is represented by the following Hermitian 
matrix: 

X=f ^0 + ^3 Xi+iX2\ ^-j^-j^g-j 



Xi — 1X2 Xq — X3 

of unit determinant. Let us now define the function: 

K,{X,z) := (p,X|p,z) := (tr(e(^)X)j = {2i{z) ■ X)-^'+^p\ (117) 

where X = eX*e~^, X = X~^det{X) if X is invertible and ■ denotes the scalar product on 
M}'^. It is convenient to use the upper half-space model for H^, in which it is coordinatized 
by (t, y),t & M+, y E C, where 
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and ds^ = dt^ /t^ + t^dy^. Then this function is exphcitly given by 



One sees from the definition that K_p(X, z) = Kp{X, z) and in the following we will denote 
Kp{z,X) ■.= Kp{X,z). 

This quantity plays a prominent role in the so-called AdS / CFT correspondence of string 
theory and is referred to in that context as bulk-to-boundary propagator. Given a real 
function in the hyperboloid we can use this propagator to reconstruct the function from its 
asymptotic data at infinity. The relation f ^ cf)^ between the realization in the space of 
functions on C and on is the explicitly given by: 

= j f^KpiX,z)fiz) = {f\p,z)Kp{z,X). (120) 

where we have written this expression in terms of the coherent states. 

Starting from the expression fll20p we can compute the asymptotics of the field from 
which we recover the function on the null cone. A non-trivial computation that is most 
conveniently done using the Fourier modes f{z) = e^iP^+p^) gives 

-{l+ip) IM. M +-i^-ip) f •^^^^ 



—ip J Ti \y — zy^^ ^Pi 

Note that the operator 

intertwines the representation (0, —p) with the representation (0, p). Moreover this operator 
is the kernel corresponding to powers of the 2-dimensional Laplacian as can be easily seen 
by taking its 2 dimensional Fourier transform, namely 



'P\dM'fi^) = J^Kp{z,y)f{y). (123) 



r(l + zp) 
Accordingly we have 

d^z 



TX 



Kp{z, y)K^p{y, w) = Ko{z, w) = ti5\z - w) (124) 



In terms of the bulk-boundary propagator the asymptotics f ll2ip reads 

Kp{X{t, x), y) tSoo —-^Kp{x, y) + ^—K^{x, y) (125) 

The advantage of the representation on the hyperboloid is the fact that the intertwining 
operator between (0, p) and (0, — p) is trivial, that is | — p, X) = \p,X), thus Kp{X,z) = 
{—p,X\p,z). The scalar product between two states on the hyperboloid can be easily com- 
puted from fll25p 

KpiX, Y) = 2 [ — Kp{X, z)K.p{Y, z) = 2^^ (126) 
J n psmhr 
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where r is the hyperbohc distance between X and Y . Since it is real we have Kp{X^ Y) = 
K_p{X,Y). One can recover the bulk-to-boundary propagator from its asymptotics, which 
can be evaluated directly 

f-{l+ip) f-i^-ip) 

Kp{X, Y{t, y)) ~i_oo —Kp{X, y) + — y), (127) 

— zp ip 

from which we deduce using (11211) 

K_p{X,z) = j '^Kp{X,uj)K.p{w,z). (128) 
The asymptotics and correspondence between (0,p) and (0, — p) therefore reads 

|p,X(t,x)) tJToo ^|p,t^(x)) + — I -p,t^(x)), 

—ip ip 
— K_p{z,uj)\p,uj). (129) 

TT 



B. Coherent states. 



Coherent states for the simple representations (0, p) are easiest to describe in the realiza- 
tion in the space of functions on the hyperbolic space H^. Such a state is parametrized by 
a point 2 G C and is given simply by: 

0|,,,)(X) ■.= Kp{X,z). (130) 

Thus, the coherent states are essentially the states |p, z) introduced above. From this defi- 
nition it is now possible to compute the scalar product between different coherent states 

{p\z\p,u) = y^^_</,|^,,>(X)0|,,,)(X) = j^^-^^K.p,{X,x)Kp{X,u;) (131) 
One can first easily compute this convolution for general points on 

[ 7|^^^-P'(>^,^)W^) = ^^^4^^p(^,^) + ^^^^W^)- (132) 
Jh» (27r)^ p^ p^ 

Using the asymptotics of the previous section we can now compute the scalar product 
between two coherent states 

(p', z\p, u) = ^SI_Ak,{z, w) + ^^P'y^ Kp{z, w) (133) 

which is the result we were looking for. 

The following identity is not hard to prove: 

/• + 00 

/ p^dpKpiX,Y) = {27rf 6{X,Y), (134) 

^0 
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where the 5- function on the right-hand-side is that on H^. Thus, we get the following 
decomposition of the identity 

lH'i = - j pMp J —\p,z){p,z\ (135) 

Note that the unity operator on is an operator that commutes with the action of SU{2), 
uIhsu"^ = l_H-3 for u G SU{2). It can be represented as an averaging operator acting on 
functions on SL{2,C) 

: L\G) ^ L\SL{2X)/SU{2)) = L^{H'), lH^F{g) = [ duF{gu). (136) 

JSU{2) 

In order to make the comparison with the Euclidean case, it is convenient to work in term of 
the representation on S'^ = SU{2)/U{1). The following decomposition of the identity holds: 



— \p,z){p,z\ = [ \p,n{z)){p,n{z)\= [ dn\p,n){p,n\ (137) 



These states are the Lorentzian analogs of the Euclidean states we used to construct the 
gravitational spin foam model. Using these states it is possible to get both the Barrett-Crane 
model as well as a new Lorentzian model rather quickly. 



C. Lorentzian gravitational model. 

Let us see how one can obtain the vertex amplitude of the Barrett-Crane model using the 
coherent states. As in the Euclidean case, the idea is to insert into the partition function of 
the BP theory the identity operator on that can be expanded in terms of a complete set 
of simple coherent states. This prescription projects onto the set of simple representations 
and identifies the simple bivectors labelled by p, n associated to neighboring tetrahedra. 
After this is done, the amplitude for a triangulation separates into pieces - the vertex and 
edge amplitudes. The vertex amplitude is parametrized by 5 x 4 quantities Ztf {ritf in the 
Euclidean case). The amplitude is given by: 

Lnew[{pf}dntf}]= / Y[dhtY[{pf,nt(^f)f\{ht^f))-^ht:(f)\pf,nt'^f)f), (138) 

iSL(2,C) J 

where as before t{f),t'{f) are the two tetrahedra that share the face / and we have denoted 
{p,z\p,w) = Kp{z,u!). Thus, the prescription for the gravitational model is to take these 
vertex amplitudes (11381) . multiply them (possibly using additional edge amplitude, see be- 
low), and then integrate over the variables describing the normals keeping these variables in 
the neighbouring simplices the same. All this is in exact analogy with what happens in the 
Euclidean case. 

Now, as in the Euclidean case, the Barrett-Crane model is obtained by integrating over 
the "normals" ritf independently, forgetting about the fact that there is a neighbouring 
simplex that one should solder the amplitude to. One obtain a model in which we have 
the characteristic Barrett-Crane decoupling of the vertex amplitudes. All this is in exact 
analogy with what happens in the Euclidean case. 
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The effect of integrations over the variables ritf is to produce the state: 

which is the SU(2) invariant state in the representation (0,p). The function (p, X|p, 0) in 
is known as the spherical function. After the integration, the vertex amplitude of the 
Barrett- Crane model becomes: 



LBc[{pf}] = [ n^^*nw'Oi(^*(/))"'^*'(/)i/'/'0)- (140) 

iSL(2,C) t J 

The quantity 

{p,0\h-^h'\p,0) (141) 

is actually a function of two points on H3 (in view of SU(2) invariance of the states |p, 0)). 
This function is the bulk-to-bulk propagator, which has been constructed in the previous 
subsection. The integrals in (11401) then reduce to a multiple integral over the hyperbolic space 
of a product of bulk-to-bulk propagators - one for each face of the 4-simplex in question. 
The quantity (I140p is thus the well-known vertex amplitude that was first introduced in [lo| . 

To finish the description of the model it remains to give the set of edge amplitudes. One 
possibility is to set these amplitudes to be just the (5-functions, similar to what we have 
done in the Euclidean context, see (jS]). Another possibility is to consider a non-trivial edge 
amplitude constructed in exactly the same way as the vertex amplitude (this is also possible 
in the Euclidean context, but we refrained from discussing it above). This edge amplitude 
is more natural from the perspective of group field theory, and is arguably the correct way 
to allow for a non-trivial holonomy between the 4-simplices. The amplitude is parametrized 
by 2 X 4 variables Zf,z'j and is given by: 

LLw [Pf. , ' 4J = / n ^Pf^ ' ^f' I ^l/'/. ' ) • (142) 

>/SL(2,C) 

The expressions (11381) . (I142p are formal, as all the above integrals over the Lorentz group 
are likely to diverge. Thus, they must be given sense via an appropriate regularization 
procedure, as is also the case for the Lorentzian Barrett-Crane model, see 1^ . It is possible 



that the same regularization procedure [10| may work for the new model. Also the amplitude 
for the whole triangulation is obtained by multiplying the vertex and edge amplitudes, and 
then integrating over the variables Ztf using as the measure the usual d'^z^f measure on the 
complex plane. This operation may again result in divergences that have to be taken care of. 
After this is done, one has to integrate over the representation labels p/, which is likely to 
result in additional divergences. All this has to be analysed in detail before it can be claimed 
that the model introduced above is well-defined. We leave this work to future publications. 

It would be important to study in more details the derivation sketched in section IVlBl in 
which the simplicity constraints are imposed by an integral over an appropriate subgroup 
of the Lorentz group. Such a generalization will give a clearer geometrical perspective on 
the constructions we have introduced. It would also be quite interesting to find a model 
that would give the Lorentzian counterpart of the Euclidean model introduced by (3]. In 
analogy to what happens in the Euclidean case it is to be expected that this other model 
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will correspond to the topological sector. It is clear that this model should be constructed 
using some other type of coherent states. Such a construction would be quite important 
for it would give a much better understanding in which sense the model introduced here 
corresponds to the gravitational sector of the theory. Again, we leave an analysis of these 
and related issues to future publications. 

Finally, let us note that unlike in the Euclidean case, there is not one, but two very distinct 
Lorentzian Barrett-Crane models. Indeed, in its most simple form the Barrett-Crane model 
is obtained by integrating over the diagonal SU(2) inserted in every strand forming the 
dual edge. This selects from the sum over representations only the simple representations 
of the type (0,p). It has been argued in STj that one can introduce the simple discrete 



representation (s, 0) by integrating over the subgroup S'L(2,R). Such a prescription can be 
adapted to our model by inserting instead of 1^3 the identity operator l^gs on de Sitter 
space. A detailed geometrical analysis of this other Lorentzian model is required. 



VII. GROUP FIELD THEORY 

The group field theory (GFT) version of the new models can be obtained quite easily, 
as our presentation above was very much in the GFT spirit. Compared to the original 
group field theory setup for the Barret-Crane model the modifications are minor in that one 
only needs to make some obvious modifications of the kinetic and interaction terms. First, 
consider the following functional on SO {4): 

T{g) = Y,^2,try,md9T,) (143) 
i 

where g = {g~,g~^) and d2jlj is the projector 

Tj := j dn |j, n) ® |j, n) (j, n\ ® (j, n\ (144) 
that we encountered above. We define the propagator to be 

K{girg^)= I dudvT\T{ugCgi\~^). (145) 



The topological spin foam model can then be obtained as the Feynman graph expansion 
of a group field theory with the following quadratic term in the action: 

c|>{gi)K{g,r9^)4>{9^)■ (146) 

50(4) 



The interaction term is given by an expression identical to the one usually used [111 . |13 . 

The gravitational spin foam model is constructed in the same way but starting with the 
intertwiner 



Gj:= j dn I J, n) ® |j, n) (j, n| ® (j, n| (147) 
instead of Tj. The Lorentzian model group field theory is constructed similarly. 
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VIII. ALTERNATIVE WAY TO IMPOSE THE CONSTRAINTS 

In the previous sections we have heavily used the coherent state techniques in order 
to constrain the discrete analog of the B field in a way that solders the geometry of the 
neighboring tetrahedra. In this section we will shown that it is possible to construct the 
same models from a slightly different, more geometrical point of view. The idea is to restrict 
the integration over the connection in a geometrical manner. 

Given a four simplex a, a discrete SO (4) connection on a is a collection of group elements 
go-j and gjj, where go-t represents the holonomy of the spin connection along the segment 
going from the center of cr to the center of a tetrahedron t, and g^j represents the holonomy 
from the center of t to the center of a face /. Our conventions are such that g^rf = g^"^, 
glj = (gj^)~^. Given a simplicial complex A, a discrete 5*0(4) connection on A is the 
choice of a discrete connection g^i, g^j for each 4-simplex. A discrete 2-form field on A is 
an assignement of a Lie algebra element X/ to each face / of A. 

Given such a discrete SO (4) connection and two form field we can write the discrete BF 
action as a pairing between the variables Xf and the holonomy (curvature) associated with 
each pair (a,/) (such a pair is called a "wedge" in the litterature). This action is exactly 
analogous to Q, except that now the sum is taken over the wedges, not (dual) faces: 

5 = 5^ tr (X;G^) , G-f ^ g'ftgtag.t'Sh^ (148) 

where t, t' are the two tetrahedra of a sharing the face /. The quantisation of this theory is 
given by the following path integral: 

The BF theory is obtained by restricting the discrete connection to be such that g^j = gj[ 
for neighboring 4-simplices. This way of writing the BF theory amplitude is equivalent to 
what we have done in the earlier sections, where we have assigned only a single quantity 
Xf (and hence Uf) per face. The main result of this section is that the new gravitational 
models considered above can be obtained from this integral by imposing constraints on the 
discrete connection alone. 

In order to describe this restriction on the connection to be integrated over we need 
to introduction an extra structure on A that we call "geometrical data". By definition 
"geometrical data" on A is an assignment of a set (e^, X/, u^) of the orthonormal frame 
at the center of a, a set of bivectors Xj at the center of each face / and a set of unit vectors 

at the center of each tetrahedron t inside a. 

Given an orthormal frame e^,A = 0, 1,2,3 it is convenient to think of unit vectors 
as elements of SU(2) ~ 5*^. The normalization condition is then tr(e'^(e'^)~^) = 26^^. The 
frame should be thought of as "sitting" at the center of a. Now, given a frame at a and a 
discrete connection in a , one can parallel transport this frame to the centers of tetrahedra 
t of cr as well as the centers of faces / using the discrete connection. Let us introduce a 
special notation for this transported frame: 

ef = (gt.e,)^ = gle^g^t, ef = (g/^e*)^ = gpfg-^, (150) 

where g"^ are SU(2) group elements representing an SO (4) one. 
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Let us now state a condition of compatibility between the SO (4) discrete connection and 
the geometrical data. Thus, recall that given a bivector X it is always possible to chose an 
orthonormal oriented frame such that X assumes the following simple form (the normal 
form) X = a A + /3 e° A where a, P are subject to the condition a > but otherwise 
arbitrary. Now given a frame and a discrete SO (4) connection we can construct from it 
the "reference" geometrical data by chosing for each tetrahedron the vector and a family 
of bivectors X/(a, /?) := a/ ej- Aej + /5/ e° Ae^. We say that a discrete connection {gat, gtf) 
compatible with the geometrical data (e^, Xj, u^) if the following conditions are satisfied: (i) 
u1 = ^ gta^a'^ (ii) gtf ■ X/ is equal to Xj(/3, a) up to a rotation fixing Ut and for a choice 
of parameters {a,f3) satisfying a > Note that these conditions are not constraints on 
bivectors Xf, but those on the frame fields and hence the connection. In particular, we do 
not assume that the bivector fields are simple. These conditions are quite natural. Indeed, 
the first one ensures that the chosen normal vector is transported consistently inside the 
4-simplex. The second one ensures that the parallel transport inside tetrahedra preserves 
this normal vector. Note that the parallel transported face vector gfj ■ Xj is exactly the face 
bivector denoted Xjj in the previous sections. The additional restriction on the parameters 
{a, (3) resolves a remaining discrete ambiguity. We can chose a dual condition /3 > |a| which 
will result in a different model. 

Our main claim in this section is that the new gravity model discussed earlier in this paper 
(corresponding to the gravitational sector of the theory) can be obtained from the integral 
(11491) by restricting the discrete connection to be compatible (in the sense just defined) with 
the geometrical data (e;^,Xy,M^) and integrating out u". Because of gauge invariance, after 
such a restriction and integration the model is independent of a choice of the frame e^. 

Proof: For definitness we fix the frame to be the canonical one e° = l,e^ = a\ 
The first condition the reads e° = g'^^Ut{g~t)~^ , and is solved by g^-t = {h^t, ha-tUt), where 
{hta, hta) is an SU(2) element that fixes the vector e° = 1. To solve the second condition let 
us denote by n/ = {n~j,nj) the S0(4) rotation which brings X/ into its normal form with 

respect to ef: n^^ ■ Xj = a A + /9e^ A with a > Applying e = (l,e) we can 
exchange the role of a and (3, thus Xf = n^e^^ ■ Xf{f],a). The second condition is then 
solved by {g^jrU'f, gt/nj) = {htf, htf){h^^j, h^^^^.){l, e) where hfj is an arbitary SU(2) element 
and 04/6 [0,27r]. 

This shows that an S0(4) connection compatible with the geometrical data {ut,Xf) is 
parametrized by two SU(2) group elements hta, htf and one U(l) element The matrix 
element entering the integral (11491) is then given by 

(j^,n^|G^|j^,n^) = (j;,j;|/i^,/i,>.,,/i?^|j;,j;)e^(^-^-^-")(<^*'/-<^*/) (151) 
X 07, -jj\h'}tK)-'hlha,u'[,ht,f\jJ, -jj). (152) 

The gravitational model is then recovered by integrating over all these data with the condi- 
tion that hi J- = hij, which implies that X^j = X^jr, where g^j ■ Xj-. Note that the condition 
^tf ~ 'Pt'f required by the geometrical condition X^j = Xfj. The integral over these 

U(l) group elements therefore imposes the simplicity condition = j~ , without having to 
put it by hand. It is relevant to note that this model does not require that u" = . If one 
relaxes the condition X^^ = X^^ one obtains the Barrett-Crane model, as we have seen in 
section HVl 

It is interesting to note that the above analysis suggests that one can impose an even 
more stringent geometrical constraint: = , which insures that not only the internal but 
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also the external geometry (the normals to the tetrahedra) agree. This gives an even more 
constrained model that would be interesting to study. However, from the canonical point 
of view the external and internal geometries are "canonically conjugate", so it is probably 
not very natural to impose constraints on both data. However, more analysis is necessary 
before a final conclusion on this can be reached. Note also that this fully constrained version 
amounts to relaxing the way the closure constraints are implemented and imposing them at 
the operator level instead of the level of the states. Indeed, for each tetrahedron there is a 
double integration f dutdu'fUtIu~[\ where J is a certain operator. This double integral allows 
to split / as a product of intertwinners, which means that the closure constraints are imposed 
strongly. If one restricts Ut = one gets only one integral left, namely / dutUtlu^^, and 
the closure constraint is imposed weakly. Note finally that we could also decide to constrain 
Ut and not Xf, the physical meaning of which would be interesting to unravel. 

To summarise, there are three different ways to impose the constraints on the geometrical 
data: (i) we can fully constrain the data by imposing = and X^j = X^j-; (ii) partially 

constrain the data by imposing either = uf or X^j = Xf^ which is the way of imposing 
the constraints we have studied in this paper; (iii) impose no constraints, which leads to the 
Barrett-Crane model. 

IX. CONCLUSIONS 

Let us summarize the results of this paper. Using Plebanski (or BF) formulation of 
gravity as the starting point we have described how the simplicity constraints that convert 
the topological BF theory into a gravitational one can be imposed using the coherent state 
techniques. We have seen how the Barrett-Crane model [5| gives a quantization of the correct 
gravitational sector of the theory but decouples the neighbouring simplex amplitudes in an 
unnatural way. We have seen how this decoupling or ultra-locality problem is resolved in 
the model of Engle et al [3] and how this model corresponds to the topological sector of 
the theory instead of the gravitational one. Finally, we have constructed a new model, 
both in the Euclidean and in the Lorentzian settings. Both models do not suffer from the 
ulta-locality problem characteristic of the Barrett-Crane models. The Euclidean model we 
obtained can be seen to correspond to the "correct" gravitational sector of the theory. In the 
case of the Lorentzian model further analysis is necessary to establish whether it describes 
the gravitational or the topological sector. However, as the Lorentzian model we constructed 
is the "delocalized" Barrett-Crane model, this strongly suggests that the model corresponds 
to the gravitational sector. 

Importantly, we have described how the Immirzi parameter 7, which at the level of the 
continuous formulation of the theory appears as a parameter in front of the "topological" 
term in the action can be incorporated into the spin foam framework. Importantly, the cases 
7 > 1 and 7 < 1 correspond to models of a rather different nature. We have seen why it was 
impossible to incorporate 7 into the framework of spin foams via the Barrett-Crane way of 
imposing the constraints. We have also seen that in the case of Riemannian signature models 
the Immirzi parameter turns out to be quantized, see (^11) and (jHS]) for the corresponding 
expressions. 

In spite of a number of interesting advances made in the present work, a great deal 
remains to be done. Thus, it is necessary to repeat the analysis of the semi-classical limit 
of the vertex amplitudes for the new models and see whether the amplitudes are dominated 
by non-degenerate geometrical configurations. It is necessary to develop the Lorentzian 
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signature model(s) in more detail, and, in particular, identify the model that corresponds 
to the topological sector. After this is done it is important to study the relation with the 
picture of quantum geometry that arises from loop quantum gravity. It is not unlikely that 
the projected spin networks of 3^ will be an important part of this relation, see also a more 



recent work [27| for some steps in this direction. It is also interesting to develop in more 
details the group field theory perspective on the new models. 

At a higher level of complexity, it is essential to continue with the program of the path 
integral quantization. Thus, in general in the path integral one has to impose not only the 
simplicity constraints, but also other second-class constraints, as well as gauge-fix the first 
class constraints (even though there may be some simplifications if one is only interested in 
the partition function, the case we have explored in the present work). So far, all these other 
steps have been missing from the framework of spin foam models of four-dimensional general 
relativity. It is unlikely that further significant progress can be made without tackling these 
issues. Another important point that has so far been missing from the spin foam framework 
is understanding in some explicit fashion how to discretize the Lagrange multipliers that 
impose the simplicity constraints. Indeed, recall that (some of) the Lagrange multipliers 
in question receive (on-shell) the important interpretation of being the Weyl part of the 
Riemann curvature tensor. This is exactly the part of the curvature that describes the 
local degrees of freedom of the gravitational field. Hence, it would be very important to see 
explicitly how the spin foam models of the type described in this work encode this Weyl 
part of the curvature. 

In spite of the spin foam quantization programme being far from complete, some prelim- 



inary physical computations of the type [38|] are possible. It remains to be seen what the 



physics described by the new model(s) is, but the absence of ultra-locality, which seemed to 
be at the root of most of the problems suffered by Barrett-Crane model(s), suggests that a 
new exciting period of the development of the subject of spin foam models is opening. 
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Appendix: Some necessary formulae 

The following two standard formulae are often used in the present paper: 

/ dgtU9)CM) = ^5™'"'<5""'. (153) 

ii+j2 



i=lii-j2| 

The last formula, in particular, fixes our choice of the normalization of the Clebsch-Gordan 
coefficients. 
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Appendix: Clebsch-Gordan coefficients 

The general formula for the Clebsch-Gordan coefficient of SU(2) (in the normalization 
that the 6'-symbol is equal to one) is given by: 

rhi.i ^ ( y,^, [r,j]{h + i2-iy. ^ i-iy {h+j + ry.{k + i-j-r)\ 

J'^™ ^ ' A(r) 4^ r\ {l-m-r)\{li-j-r)\{l2-l+j+r)V 

(155) 

where the sum is taken over all values of r such that the arguments of the factorials are 
non- negative integers (zero is fine), also m = j + k, and 

Ail) : = ((/i + I2- /)!(/ + /i - kV-il + I2- hV-ih + I2 + I + iV-f^ . 
For purposes of this paper we need to specialize to h = h = L. The formula becomes: 

(2L - 1)<(L - mL - k)\ 



(-1)' (L+j + r)!(£ + i-j-r)! 

r! {I — m — r)\{L — j — r)\{L — I + j + r)\ 



The special case of this formula that we need for purposes of this paper is that of one of the 
vectors being the highest weight. Thus, we choose j = L. It is easy to see that in this case 
the requirement that the argument of the factorial {L — j—r)\ is non- negative truncates the 
sum over r to the single term r = 0. The formula becomes: 

^LLi _( {2L)\{L ~ m + L + k)\ y/' 

LkL+k \{2L-l)\{2L + l + l)\{L + k)\{l-L-k)\) ■ ^ ^ 

Some special cases of this formula are: 

can = (159) 

rLLi _ ( m\ (2L)! 

" [i2L-l)\i2L + l + l)\J ■ ^^^"^ 

It is interesting to compute the asymptotics of the last coefficient in the regime 2L — 1 » 1. 
Using the Stirling formula one obtains: 

ci:'J,o^ , . ^1 . . (1611 



V2L + / + 1 \l + ^J (l-(^)' 



Now, in the limit where L — 00, 1/2L = x fixed and strictly smaller than one we have an 
exponential decay 

CL"« - 4^ ,. 1. . f^^^) " . (162) 



2L(1- 1)1(1 + i)J V(l + a:)' 
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since one can easily check that the term inside the exponent is smaller than 1 when x < 1. 
When X is itself infinitesimal we have 



(l + a;)i+^ 



_ .2 



e"" . (163) 



The general formula fllSSp also gives the set of coefficients relevant for the non-trivial 
Immirzi parameter case. These are given by: 

^.-r. ( (2J")! (2r)! 
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